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Abstract 


The world has seen an increase in the popularity of renewable-based 
energy systems. However, because of their erratic nature, fossil fuels 
cannot be entirely phased out. Utilizing carbon dioxide in greenhouses close 
to power plants and generating extra power electricity from dissipated 
thermal energy are appealing strategies that cause us to have a cleaner 
environment and affordable power electricity. To exploit the high enthalpy 
of expelled gases from a gas turbine, the gas turbine was coupled with an 
absorption refrigerant cycle and an organic Rankine cycle and analyzed. 
Additionally, a deep artificial neural network method was used to calculate 
all functions in a rapid optimization and accurate 4E analysis. The carbon 
dioxide emissions subsided from 0.9183 to 0.41 kg CO2/s, and the 
optimization time improved 257 times after using the trained machine 
learning model. An adsorption technique for CO2 separation is used for the 
proposed system because of its lower energy consumption. The maximum 
exergy and energy efficiencies were attained at 36.71% and 49.2%, 
respectively, and parametric studies are being assessed. Due to increases in 
harvests and efficiencies, the net annual interest has increased by 21.8%, 
up to 22.5 million dollars. 


Keywords: Toward Zero Carbon; Environmental impacts; Power and air- 
conditioning system; Greenhouse; Deep learning; Absorption Refrigerant 


cycle; 


1. Introduction 


How to reduce air pollutants is one of the most important factors in 
energy systems. Most academics try to use techniques to absorb pollutants 
like carbon dioxide and methane generated in this setting [1]. Taking 
advantage of the carbon dioxide generated in greenhouse environments for 
plant photosynthesis is one of the most common strategies [2]. It is possible 
to assess a self-demanding combined heat and power system coupled with a 
greenhouse to absorb carbon dioxide [3]. Utilizing hot flue gas for heating, 
cooling, or dehumidification while simultaneously generating electricity on- 
site, combined cooling, heating, and power systems can greatly increase 
energy efficiency [4]. According to the effects of renewable energy, natural 
gas is the most environmentally friendly fossil fuel [5]. Utilizing hybrid fuel 
cell systems with gas turbines can increase energy generation while 
reducing polluting emissions by up to 47%. To calculate how much the 
cooling, heating and power system in combination with a greenhouse may 
reduce the carbon dioxide footprint, it is possible to compare conventional 
and organic harvests with the CHP [6]. CHP affects both conventional and 
organic crops equally. Nearly all energy systems emit carbon dioxide; 
hence, researching several multigenerational systems with diverse 
configurations is essential in terms of energy, efficiency, and economics. In 
a real-world case study, these three assessments are carried out 
independently and concurrently for three different types of gas engines, 
gasoline engines, and gas turbines as power units [7]. Greenhouse crops are 
often grown using the heat and carbon dioxide that energy systems create 
[8]. Thus, it is necessary to determine the needed carbon dioxide and 
heating load. Additionally, the greenhouse's needed cooling load during the 
summer might be met by using the heating load to power an absorption 


refrigerant. Consequently, an absorption refrigerant simulation [9,10]. 


The energy and exergy studies of each component in such systems are 


strongly advised to shed additional light on the subject. Calculating the heat 
transfer of heat exchangers between the energy production unit and the 
absorption refrigerant cycle or between the gas engine and the greenhouse 
atmosphere significantly impacts the model's development [11]. As a result, 
a portion of the power produced by the Bryton cycle or other power units is 
utilized to power pumps and compressors as well as to illuminate the 
greenhouse [12]. Consequently, both the lightning process and other system 
components should be represented [13,14]. Additionally, the greenhouse 
environment's heating value and cooling load should be determined 
appropriately depending on the appropriate temperature for crops growing 
there. In this context, it is also possible to determine the amount of energy 
required to compensate for temperature changes in warm or cold seasons 
[15]. In such systems, irrigation is also taken into account, and the water 
requirements for crops grown in greenhouses should be determined [16]. To 
gain a thorough understanding of the subject, a practical, real-world case 
for modeling of irrigation and determining the amount of water required for 


specific types of harvests was employed [17,18]. 


The energy industry has noticed the appeal of carbon dioxide capture, 
which has been studied from many angles, including energy, efficiency, 
economy, the environment, and even the social impact on urban dwellers. 
Although there are several techniques to absorb carbon dioxide released 
from power plants' energy production units, examining these strategies 
from an economic perspective would be preferable. Therefore, it might be 
helpful to understand these processes to estimate the effects of coupling a 
Brayton cycle with an absorption cycle using a Mono-Ethanol-Amine (MEA) 
solvent [19]. It is also modeled and examined to integrate gas turbines with 
molten carbonate fuel cells in power plants to lessen or eliminate carbon 
dioxide produced by hybrid systems [20]. In such a system, energy and cost 
analyses are done while simulating a molten carbonate fuel cell. 


Furthermore, the Biomass Integrated Gasification Combined Cycle 


technology allows for the collection or regulation of carbon dioxide 
emissions [21], and economic assessment might be used to gain a thorough 
knowledge of such technology to be taken into consideration as a different 
way to capture carbon dioxide. The exploration of more energy- and cost- 
efficient arrangements and procedures might be useful because these 
multigenerational systems release a significant quantity of carbon dioxide 
[22]. In this context, an integrated system with a carbon dioxide extract and 
storage unit, an Organic Rankine Cycle, and a chiller cycle may be 
evaluated from the standpoints of energy, efficiency, the economy, and the 
environment [23]. Additionally, phycocapture, which may be considered a 
carbon sink in metropolitan areas, is a unique way to reduce or capture 
carbon dioxide because greenhouse gases like carbon dioxide might 
profoundly affect many aspects of human existence [24]. By using 
microorganisms like photosynthetic bacteria and microalgae, carbon dioxide 
will be trapped; for these organisms, carbon dioxide serves as a source of 


carbon and aids in their growth [25]. 


When designing a thermodynamic and economic model, one approach 
is to use a MATLAB function for analysis. However, maximizing the 
objective functions can be challenging due to the complexity of the model 
and the need to balance multiple factors such as efficiency, cost, and 
environmental impact. An Artificial Neural Network (ANN) or deep 
learning-based methods may be useful in particular situations since the 
complexity of the task can be increased by the quantity of objective 
functions and decision variables. Machine learning-based issues may be 
solved using a variety of techniques and algorithms [26], and any of which, 
depending on the situation, the decision-making factors, and the 
accessibility of the data collection, might be utilized for a unique case study 
[27]. Additionally, comparing metaheuristic algorithms with machine 
learning can help us comprehend machine learning issues in depth [28]. 


These algorithms have been contrasted with classical optimization 


techniques like the Archimedes optimization method [29]. S. Ikeda and T. 
Nagai have introduced a unique optimization approach since it is hard to 
compute and go forward using traditional algorithms when the number of 
inputs is vast [30] to cut down on the input variables. Additionally, 
metaheuristic algorithms can be integrated with a multi-layer reinforcement 


learning technique [31]. 


The vast majority systems of multi-generational have not stated a 
thorough examination; among them emphasize the thermodynamic 
viewpoint versus others concentrate on the economic field. Here, we 
discovered that all thermodynamic viewpoints, economic theories, and 
environmental theories require a thorough investigation of CCHP systems. 
This paper has presented a study and optimization of CCHP and greenhouse 
units to utilize the CO2 in a power plant's exhaust gas. The system includes 
a Brayton cycle with an intercooler exchanger for producing electricity and 
extracting heating energy, an Absorption Refrigerant Cycle for cooling or 
heating a greenhouse unit, an Organic Rankine Cycle for producing 
additional electricity, and a power management unit. In order to increase 
the performance of the integrated system it is tried to formulate a multi- 
objective optimization problem and came up with the optimum decision 
variables giving the Pareto-front of the multi-objective function optimization 
problem. To reduce the computational demand and runtime of the problem, 
a Deep Neural Network is used, which is trained based on the results of the 
model that is derived for the integrated system. It solved the optimization 
problem through a detailed analysis of various designs based on which the 
CO2 capturing system can operate and showed the benefits of the ultimate 
design based on the 4E analysis of the system. Also, the study did not focus 
on separation techniques or catalytic converters, but on separation 
techniques; the adsorption process is preferred for this system because of 
its lower energy consumption and financial costs. In addition, one of the 
most important objectives of the presented study is CO2 capture. Thus, 
despite the fact that the ORC would need a higher temperature than the 
absorption chiller, providing cooling and heating energies and other 
conditions for the greenhouse are priorities to capture CO2. For the 
greenhouse to always produce crops at optimum conditions and consume 
CO2 to clear the environment, it should provide various types of energy, air 
conditioning, and all other requirements for various conditions. For these 
reasons, it was decided to adjust the priority for the greenhouse. 


We looked at a genuine power plant and a greenhouse that is either 
already there or will be constructed soon to lessen the environmental 
effects of the power plant by absorbing the carbon dioxide it produces. In 
this case, we examined economic and environmental factors when modeling 
the complete power plant and CCHP system in terms of thermodynamics. 
Table 1 lists the areas the present study fills up where other earlier 


investigations left off. 


Table 1: Identifying gaps in the literature and contrasting them 
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2. System description 


The system, illustrated in Figure 1 consists of a Brayton cycle, an 
Organic Rankine cycle, a greenhouse, an air conditioning system and chiller 
cycle, and will be evaluated as part of the project. The project's main goal is 
to employ a greenhouse to capture CO2 produced in power plants in order 
to reduce their impact. The photosynthetic process is carried out by 
greenhouse plants using carbon dioxide and sunshine. So, by supplying the 
necessary standard carbon dioxide, the suggested method has the potential 
to lower air pollution and enhance greenhouse products. Watermelon and 


tomatoes are two of the plants planted in the greenhouse. 
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Figure 1: A diagram of the suggested system 


The detailed system functioning is shown in Figure 1. Air enters point 1's 
initial compressor, which then compresses the air. To reduce the power 
consumption of the compressor from points 2 to 3, the temperature of air is 
decreased by a heat exchanger, which for the cooling process it uses cold 
water to heat up for domestic water applications. The second compressor is 
then filled with the cooled and compressed air. Fuel and compressed air 
enter the nonadiabatic combustion chamber with the efficiency of 90% at 
points 4 and 5, respectively. The combustion process occurs in the 
combustion chamber, and the high-energy flue gases exit the combustion 


chamber at point 6 and enter the gas turbine to provide the needed power. 


Before transferring the electricity to the grid, the system deducts the 
necessary power for the fluid delivers. The exited gases still have a 
significant quantity of energy at point 7. After some related chemical 
processes in the urea tank and catalytic converter to remove toxic species, 
the flue gas reaches point 9 to transfer the energy to heat exchanger of 
absorption refrigerant cycle. A catalytic converter is a device that is used to 
reduce the emissions of harmful pollutants from the flue gas of a 
combustion engine. In a gas turbine cycle, the catalytic converter can have 
a positive effect on the efficiency of the cycle by reducing the amount of 
pollutants that are released into the ambient and the combustion again. The 
present study assesses an oxidation catalytic convertor to reduce pollutants 
and practical issues. The process of eliminating harmful pollutants is 
achieved through a chemical reaction that performs on a catalyst surface 


within the converter. 


Exhaust gases exit the heat exchanger at point 10. Whether or not point 
10 has enough energy to run another subsystem depends on the time of 
year and the gas turbine's power output. Fortunately, point 10 frequently 
has enough energy to take into account using an Organic Rankine cycle to 
produce more electricity. As a result, the expelled exhaust at point 10 
transfers to a different heat exchanger, where it heats the Organic Rankine 
cycle fluid before leaving at point 11. Carbon dioxide is provided to the 
greenhouse at point 13 using carbon dioxide separation procedures, while 
other subsiding flue gases exit the system at point 12 through similar 


operations. 


The CO2 separation process at the end of a gas turbine cycle involves 
several steps, including capturing the exhaust gas through a flue gas duct 
from the turbine, cooling to remove moisture and other impurities, 
separating the CO2 from other gases using various separation techniques 


such as absorption, adsorption, or membrane separation, compression to 


increase its density for transportation, transportation to the greenhouse. 
There are different ways to separate carbon dioxide from the flue. The gases 
resulting from combustion have a high temperature, their pressure is low, 
and the concentration of carbon dioxide is low. For this reason, absorption 
and adsorption methods are suitable for separating carbon dioxide for post- 
combustion processes [35]. The adsorption process can lower the energy 
and financial costs associated with CO2 separation from the flue. Generally, 
the CO2 adsorbent has to be highly selective, have a large adsorption 
capacity, have acceptable adsorption or desorption kinetics, and have 


strong thermal and mechanical stability. 


The air conditioning subsystem encounters two basic scenarios; in the 
winter, the greenhouse requires heat. The greenhouse receives the heat 
that was transmitted from the first heat exchanger. Another circumstance is 
when the greenhouse needs a cooling load during the warmer months. The 
first condition occurs at points 14, 15, 16, and 17. In warmer months, the 
transferred heat drives an absorption refrigerant cycle from the first heat 
exchanger. For point 22, the steam generator receives this heating load, 
which leads the refrigerant fluid to become a saturated refrigerant. At point 
20, the weak solution leaves the steam generator and has a lot of energy 
that might be utilized to boost the energy of the strong solution. Thus, it 
travels to point 26 after going via a solution heat exchanger. Then exits a 
liquid expansion valve, and until it reaches point 27, its enthalpy and 
pressure are decreasing. The refrigerant enters point 23 as a saturated 
liquid after passing through a heat exchanger to condense and descend its 
energy. After that, it moves through another expansion valve to arrive at 
point 24 in a two-phase condition and at a reduced pressure. This two-phase 
flow is prepared to draw energy from the surrounding environment, namely 
the greenhouse. In this case, the refrigerant travels from an evaporator to 
saturated vapor, where it delivers the greenhouse's necessary cooling load 


before reaching point 25. Both the refrigerant and the weak solution enter 


the absorber at point 28, where they create the strong solution and release 
heat into the environment. At point 28, a solution pump presses the 
saturated, strong solution under pressure. To reduce energy usage in the 
steam generator, the strong solution at point 29 gets heat from the solution 
heat exchanger. As a result, the condition at point 21 is getting closer to the 


ideal state of a saturated liquid. 


The second heat exchanger's heat transfer enables the organic Rankine 
cycle HTF to be heated up at point 30, where it enters the organic Rankin 
cycle turbine and exits at point 31. Through the use of a condenser, the 
ORC fluid is condensed and chilled until it reaches point 32 as a saturated 
heat transfer fluid. The ORC condenser's cooling procedure also makes use 
of residential hot water (points 38 and 39). The ORC fluid departs the pump 
at point 33 to obtain the necessary energy from the organic Rankin cycle 
heat exchanger after the ORC pump makes up to the prior pressure of the 
fluid at point 32. The heat losses have a significant impact on the 
temperature of the flue gas, because they are highly related to radiation, 
and thermal radiation is proportional to 7*. For this reason and modeling 
the lost heat, decreased effectivenesses for all heat exchangers are applied, 


which is presented in Table 5. 


3. System modeling 
Modeling of various aspects is presented in this section. 


3.1. Modeling of energy and exergy 


3.1.1. System thermodynamic modeling 


Equations for the balance of mass, energy, and exergy are expressed as 
follows [36-38]: 


» M,= ST gy, (1) 
SEn=SEn,, (2) 
SEX,— > EX = EX soot (3) 
The following equations might be created by expanding these equations: 


Q+ Ym, [Dn —h,|= W+ Stitt tev ~ hy| (4) 


T\ 3 ; 
s1-P)e,- W+ SM, €X,— S Mout @ Xour= FE Xaest (5) 
k 
Where Ex is calculated by multiplying the mass flow rate by the specific 
exergy. 
Ex=mex (6) 


The four basic categories of exergy are potential, kinetic, chemical, 
and physical exergy. Potential and kinetic energy may be insignificant in the 
current model because there are no changes to any of these quantities [39- 
42]. 


ex=ex'+ex™ (7) 


The following formulae can be used to compute the chemical and physical 
exergy [43,44]: 


ex”’= h—h,— T,|s—5| (8) 


The following might be used to calculate the particular physical 


exergy of organic streams for constant specific heat capacity: 


ex"=C,| T-T,- Tyln (9) 


L 
To 


The specific chemical exergy is also calculated by the following 
equations, which provide a unified correlation for estimating the specific 


chemical exergy of solid and liquid fuels [45]. 


€Xoy= 363.439 C+1075.633 H—- 86.308 0+ 4.14.N+190.798S—21.1A 
(10) 


Where C, H, O, N, S and A stand for carbon, hydrogen, oxygen, 
nitrogen, sulphur and ash contents of the fuel, respectively, in weight 


percent and dry basis. 


Physical and chemical exergies for gas streams might be calculated as 


follows: 

ex"=C.| T-T,—Tyn | -||+R Tn | + (11) 
sd T) Py 

oD Oy,é+RT, > Oy,ln \vi (12) 


It is important to first identify the molar fractions of different species 
that are exhausted from combustion reactants to calculate the rate of 
chemical exergy. The fuel type is a kind of organic fraction of municipal 
solid waste [46]. Assessing the combustion of an air-to-fuel mass ratio of 13, 
the molar fractions of the components of expelled gas are computed as 
follows [47]: 


0.73704 CH,+0.24176 CO,+i0.00481 O,+0.01630 H, O+5.69698 N, 
+41.57493 O,+0.00229 CO,+0.14533 H,O-0.98109 CO, 
i+1.63571 H, O+0.10566 O,+5.69698N, 


(13) 


The following equation is utilized to compute the rate of exergy rate 


because of heat transfer from or to each component of the system [48]: 


E xX; = Ons T 


1-2 (14) 


To calculate the amount of energy lost or absorbed from the 
environment, equation (14) was used. The balance equations for all systems 
are found in equations (1-14), and Table 2 contains the governed relations 


for all components. The exergy of the fuel is equal to 831200 kJ/kmol, and 


its molar mass is 235.25 kg/kmol [47]. 


Table 2: Formulae for each system component's energy and exergy 


components 


Energy analysis 


Exergy analysis 


Intercooler HEX 


Compressor I 


Air compressor II 


Combustion chamber 


GT 


Heat exchanger I 


Steam generator 


Refrigerant absorber 


Condenser of chiller 


chiller evaporator (for 
summer) 


ARC pump 


ARC expansion valve 
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The empirical correlation proposed by Macchi and Astolfi provides an 
estimate of the turbine efficiency in a Brayton cycle based on specific 


parameters. Here is the formula for turbine efficiency [49]: 


Niurp= at bx( i Tota) + cx( eee dais) + dx( Pham Proia) 


Where T;,,, is the temperature of the hot gases at the turbine inlet (in 
Kelvin), To1q is the temperature of the cold gases at the turbine outlet, P,,; is 
the pressure of the hot gases at the turbine inlet (in bar or absolute 
pressure units), P..,, is the pressure of the cold gases at the turbine outlet 


(in bar or absolute pressure units) a, b, c, and d are coefficients determined 


empirically based on the specific turbine design and operating conditions. 


The following equations are used to compute the absorption 
refrigerant cycle's coefficient of performance for energy and exergy 


efficiency [47]: 


COP=-— Quva 


Qcen + Waas 


(15) 


W wect ag 


=i ae 1 
Nen Q.+ WY cca ( 6) 
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(17) 


The combustion chamber lasts heat through radiation and convection 
mechanisms. Non-adiabatic combustion decreases efficiency due to heat 
loss from the combustor. Thus, it is not adiabatic, and modeling it as non- 


adiabatic as follows, which is also presented in Table 2: 


Qhss= O ae Qrad 


The convective heat transfer rate can be calculated using Newton's law 
of cooling: 
Qconv=A xA x( i Tint) 


Where, Q..n, is the convective heat transfer rate (in watts), h is the 
convective heat transfer coefficient, A is the surface area through which 
heat is transferred, J, is the surface temperature of the combustion 


chamber, and T,,,; is the temperature of the surrounding fluid. 


The radiative heat transfer rate can be calculated using the Stefan- 


Boltzmann law: 


Qrad=EX OX A X| T;- Ti 


Where Q,,, is the radiative heat transfer rate (in watts), € is the 
emissivity of the combustion chamber surface, o is the Stefan-Boltzmann 


constant, and A is the surface area of the combustion chamber. 


One of the best options for CO2 separation from flue gases is 


temperature swing adsorption (TSA), which uses a solid adsorbent like 
zeolite 13x and has been applied for the present study. There are some 
advantages to it, such as high CO2 adsorption capacity and selectivity under 
flue gas conditions, regeneration achieved by simple heating instead of 
pressure cycles, thermal swing providing a strong driving force for 
desorption, relatively low energy needs for the process, and a widely proven 
adsorbent with predictable performance. The CO2 adsorption slightly 
reduces efficiency due to the energy for adsorbent regeneration, but this 
process has a positive impact on environmental issues. The adsorbent is 
packed into fixed beds. While one bed adsorbs flue gas CO2, the saturated 
bed is regenerated by raising the temperature using a hot sweep gas. 
Multiple beds are cycled to provide continuous operation. TSA gives high- 
purity CO2 without needing additional displacement gases. The process is 
also relatively simple to operate. To decrease the temperature of the CO2, a 
cooling load is needed, which reduces total performance a little but has a 
significant positive effect on the environment. In this study, the cooling load 


is provided by an absorbent refrigerant cycle. 


In a gas turbine cycle or Brayton cycle, the catalytic converter is used to 
clean the exhaust gases before they are released into the environment. The 
gas cleaning process of a catalytic converter involves steps, including I) 
Oxidation: The first step in the gas cleaning process is oxidation, where any 
unburned hydrocarbons and carbon monoxide are oxidized to CO2 and 
water vapor. This is done by passing the flue gases through a catalyst bed 
that includes platinum or palladium. II) Reduction: Almost all of the 
nitrogen oxides (NOx) present in the exhaust gases are separated into 
nitrogen and oxygen. This is done by passing the exhaust gases through 
another catalyst bed that contains rhodium or platinum. III) Adsorption: 
Almost all of the remaining pollutants, such as sulfur dioxide (SO2) and 
particulate matter, are adsorbed onto a surface of activated carbon. IIII) 


Regeneration: The activated carbon or zeolite is regenerated by heating it 


to a high temperature to burn off any adsorbed pollutants. 


Oxidation catalysts are utilized to convert carbon monoxide and 
unburned hydrocarbons into carbon dioxide and water. The present study 
utilizes this kind of catalyst. The oxidation catalyst increases efficiency by 
combusting any unburnt hydrocarbons in the exhaust, recovering some 
energy. Precious metals like platinum, palladium, and others create sites 
that facilitate combustion reactions. Diesel engine exhaust systems often 
employ oxidation catalysts since diesel combustion generates amounts of 
nitrogen oxides (NOx). Here is a description of modeling an oxidation 
catalyst: 


CO+s 0; +CO, 


C,H,+(x+yl4) O, >xCO,+(y/2) H,O 


The reaction rates can be expressed as power law functions based on 


partial pressures: 
rCO= kCO-PCO” : PO; 
rHC= kHC:-PHC -PO% 


Where k is the temperature-dependent rate constant, and m, n, p, and q 


are the reaction orders determined empirically. 


The component molar flow rates are balanced between the inlet and 
outlet: 


FCO,€+rCe0O:-Wcat= FCO, out 


FHC ,¢€+rHC:-Wcat=FCO,, out+ FH, O, out 


Where Wat is the catalyst mass, and F is the molar flow rate. 
The conversion efficiency for each pollutant is: 
XCO=(FCO,¢€-FCO, out)/ FCO,€é 
XHC=(FHC,¢€-FHC, out)/ FHC,€é 


The model uses the kinetic rate expressions within the component mass 
balances to predict the CO and HC conversion across the catalyst. The rates 
depend on temperature, pressure, and reactant concentrations in the 


exhaust stream. 


3.1.2. Lighting modeling 


Light is required for all greenhouse products to be photosynthetic. 
Sunlight throughout the daytime provides this necessary illumination. It 
might not be enough. Thus, a backup light source from an artificial source, 
such as LED lighting, should be considered. Right on top of the plants, 
about two meters above them, are LED lamps that provide the necessary 
light for photosynthesis to occur at night or on overcast days. In the 
literature, modeling and estimation of the quantity and kind of lighting have 
been discussed [13,50]. The ideal quantities and configuration of the bulbs 
can also be taken into account [51-53]. A greenhouse's lighting system must 
be designed economically. Each lamp's optimum illuminating area for 
optimal plant development is noted on the packaging as Ajamp [54]. 


A cy 


NOL= (18) 


lamp 


where NOL is the number of needed lamps, Ag, is ground area of the 


greenhouse. 


In the current study, 14-Watt LED bulbs are taken into account for 
greenhouse illumination, and the energy needed for lighting the system may 


be calculated using the formulae below: 


NOL x14 
W ou tight= 1000 


(46) 


where Wz, tight is the amount of energy needed to illuminate a greenhouse 
(kW). 14 is the lights' nominal consumption power (in W), and 1000 is used 
to convert W to kW. 


3.1.3. carbon dioxide modeling 


The required carbon dioxide for typical and standard conditions is 200 
and 900 (l.l'), respectively [55]. Thus, the power plant should make up for 
the shortage of carbon dioxide in the air. The following equation yields an 
evaluation of carbon dioxide capture from gas engines or decreased carbon 


dioxide in the environment: 


Menge Acex Rue CO, (19) 
EMI 0, extra=Mc0,~ M req,co, (20) 
3.2. Modeling of economic and environmental analysis 


The 15-year system lifespan and the effective interest rate are the 
foundations of the economic analysis. The cost index of each component has 


been taken into account [56]: 


Cost index forthe reference year 
Cost index for the year when the original cost wa: 


Cost atthe reference year= Originalcost x 


(21) 


The following formulae could be used to calculate the cost rate of the 
k-th component within the system, which contains both the component's 


capital cost and operation and maintenance costs [57,58]: 


Z aL +2,” (22) 

_ Z,xCRFX® 

Z,=— G (23) 

crre ite (24) 
(1+i)"-1 


where n is the life cycle of the presented system and / is the interest rate. 


The following equations estimates the overall cost rate: 


C-= CX I ino) X 3600 (25) 
Conv= Coo, X Neo, ¥ 3600 (26) 
Cidh cost 2 OZ, + C+ Cony (27) 
Cou, ine= Cou, ine* Acu (28) 
Crtec, inc= Celec * Wit x 7500 (29) 


where M,,.) is the fuel mass flow rate and C; is the fuel unit cost. Ca. stands 
for the environmental cost rate; Cco, for the CO2 unit cost; and co, for the 
mass flow rate of CO2 emissions. The yearly revenue from selling 
greenhouse goods is Coy i, ($) while the annual revenue from selling 
electricity is Cytec, inc ($). Cox, inc ANA Cyye, stand for electricity price per kWh 


and greenhouse revenue per square meter, respectively. Table 4 displays 


the capital cost functions for each component of the system. Some input 
variables for the presented correlations in Table 3 should be defined and 
calculated, such as heat transfer rates, power rates, and areas. The heat 
transfers and power rates are calculated based on previous analysis, and 
the heat transfer area should be calculated. The design of a heat exchanger 
typically involves the calculation of the required heat transfer area based on 
the heat transfer rate, the heat transfer coefficient, and the temperature 


difference between the two fluids in the exchanger. 
OFUXA XAT ip (30) 


Where Q is the heat transfer rate, U is the overall heat transfer coefficient, 
A is the heat transfer area, and AT;yrp is the temperature difference 


between the two fluids as follows: 


Where AT, is the temperature difference between the hot fluid at the inlet 
and outlet, and 4T, is the temperature difference between the cold fluid at 


the inlet and outlet. Table 5 also provides all of the constant variables 


related to economic and environmental analyses. 


CO2 emissions contribute to climate change and its associated 
impacts and damages. In addition, CO and NOx emissions cause direct 


health impacts and environmental damage on local and regional scales. 


Table 3: Emission costs 


Cost 
Pollutant (USD/kg 
) 


From 
climate CO2 0.05 - 0.1 
change 
Health and 
environment CO 1-11 
al 
Health and 
environment NOx 5-14 
al 


Table 4: The suggested system's component capital cost functions 


Component Cost function ($) Refere 
nce 
Intercooler HEX log yo Zi= 4.3247 —|0.3030 log,o| Al|+/0.1634 (logyg| / [56] 
‘ ; z= 39.5m, || P, - P, 159] 
ir compressor ACT 0.9— Tac P, P, 
: P P 
Air compressor I AC™ oe outs = Lal 3 [59] 
0.9— Nac]\ Ps Ps 
25.6 1 
oo=| "+ |[1+exp|0.018 7,-26.4)] 
Combustion chamber 6 [5 9] 
0.995-— 
P, 
266.3 5 P. 
Heat exchangers I Lyx = 2143 <A [60] 
Heat exch 2 
“forwister) Zig 2143 (Aga! (60) 
Absorption 


Z=1144.3 x(Q 


0.67 
oa 


refrigerant fluid 


[7] 


ORC heat exchanger Zorc, Hx= 2143 x(A eg [60] 


ORC turbine Zia =4750| w,\” [61] 
ORC condenser Zorc,con =I 16.62 | Abva ler [61] 
ORC pump Zig g=200| Wel [61] 
Greenhouse Zoreenhouse =Asgreenhouse X Coreenhouse [ 6 2 ] 
Catalytic converter 5000 [ 62 ] 
Urea tank 3000 [62] 
Lamps 25 X Niamps [ 13 ] 
0.67 
Air blower Z .p=91562 x as 57 
- AB 455 [57] 


Table 5: Constant values of economic analysis [57-59,63] 


variable Value Unit 
N 7500 Operating hours 
@ 1.06 
Cco, 0.024 $/kg 
Cen 30 $/m? 
i 12 % 
n 15 Operating years 
Cruel $/kg° 
Colec 0.24 $/kWh 
Cor, inc 107.639 $/m? 
3.3. Deep neural network approach 


A deep neural network is applied for this study because, after training 
the model, the obtained precision is greater than that obtained by 
conventional neural networks, working safely with data is easier, and the 


data requires less pre- and post-processing. Thus, despite the fact that 


neural networks, but failed to produce a suitable network. Applying the 
deep neural network causes less pre-processing of data because the data 
vary over a large range. However, one of the significant limitations of deep 
neural networks is their extensive parameterization, which can result in 
over-training and limited generalization generalizability. There are several 
problems related to not scaling data for neural networks, including 
producing very high bias, which causes the model to be unstable, slow 
convergence, unstable gradients, etc [64]. Large weight values can make 
the model more sensitive to small changes in input values, which can lead to 
instability and poor performance. On the other hand, scaling the data does 
not always lead to perfect precision, and during the scaling process, 
mistakes might be made due to several sequences and details. Deep neural 
networks do not have these weaknesses. Therefore, the training is 
performed through a deep neural network. In addition, a deep belief neural 
network can learn hierarchical representations of data, which makes them 
well-suited for tasks such as feature learning and dimensionality reduction, 
which is a benefit over conventional neural networks for the mentioned 


weaknesses [65]. 


Multi-layer training of various descriptions utilizing a hierarchical low- 
level feature-based framework is referred to as deep learning [29]. The most 
notable benefit of this learning process is representation learning, which is 
a method for automatically extracting features from low-level inputs. With 
the help of Bolter Mann’s limited layer learning machine, deep belief 
networks are constructed. The Deep Belief Network, which may reveal the 
deep characteristic of the primary data, is a multi-layer of Boltzmann 
limited systems. The learning process is started by a low-level Restricted 
Boltzmann Machine (RBM), which receives inputs from the DBN and 
gradually accounts for the grading, until finally the RBM gets the highest 
level comprising DBN outputs will be learned as well. This is because the 


DBN seeks to increase the likelihood of a learning input and output data. 


This process offers a useful learning path by including a variety of 
straightforward models. The result of the first layer is then utilized as the 
input for the second layer, and so on, until all the layers have been trained 


and the outputs depart from the real output data the least. 


The following are the most crucial Deep Belief Neural Network formulas 
[29]: 


AL, Liar, xe Feet (32) 

F,= Eee (33) 

Ei Ly LiJ=-2 OBL, 2, O Bla) OL, L, wy (34) 
i= J= I,J 


Aw,=E{L,L,)-En 


i) (35) 


where L, and Ln, respectively, stand for the visible layer and the hidden 


layer. For the purpose of establishing the joint distribution function of the 
hidden unit I and the visible unit j, energy function p(L,.L,| is applied. F, is 
a partition function that can sum up pairs of units that are both hidden and 
visible. The energy of the joint configuration of units, denoted by F' | | pe | in 
Eq. (29), is derived from Eq. 30. Biases in visible and hidden units are B;, 
and B, respectively. The weight of hidden and visible units is distinguished 


by w,. The RBM weights are updated in Equation 31. 


We created a MATLAB program to provide 2000 data inputs for each 
choice variable and 2000 data outputs for each objective function in order 
to train the DNN. Despite the availability of the precise function, 
calculations take a long time to perform, especially during the optimization 


phase. Therefore, in order to shorten the computation and optimization 


processes' run times, the generated data were utilized to train the DNN. 


Figure 2 shows the DNN schematic for the current model. 
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Figure 2: Schematic of the Deep Belief Neural Network Structure 


4. Results and discussion 


It is true that certain input variables are provided for model 
computation, while others need to be approximated or assumed. The input 
variables or data are shown in the tables and graphs below, such as Figure 
3, showing Tehran's weather and temperature conditions. These data were 
averaged over a period of ten years and used in this study as inputs for the 


ambient temperature. Table 6 lists the additional necessary input variables. 


Teheran, Iran Climate Graph (Altitude: 1191 m) 
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Figure 3: Tehran's mean temp over the course of ten years, according to the World 
Weather Organization [66] 


Table 6: Necessary input variables for 4E analysis 


Variable Amount Unit 
Efficiency of Gas turbine exergy, er 90 % 
Efficiency of combustion chamber, ncc 90 % 
Low Heat Value, LHVcua 500,50 kJ/kg 
Low temperature water, Tootd, domestic 25 °C 
High temperature water, Thor, domestic 70 °C 
Gas turbine exergy efficiency, Ner 90 % 
Combustion chamber efficiency, fcc 90 % 
Compressors' exergy efficiencies, Mcomp 85 % 
Organic Rankine cycle pump exergy efficiency, BG % 
T]ORC, pump 
Organic Rankine cycle turbine exergy efficiency, BG % 
TIORC, tur 
Absorption refrigerant solution pump exergy a os 
efficiency, Marc pump 
Chiller cycle upper pressure, Pyen 1350 kPa 


Chiller cycle lower pressure, Paps 250 kPa 


HEX effectiveness, € 
Pressure drop, APcc 
HTF of absorption 
Strong solution mass fraction, X;; 
Weak solution mass fraction, Xs 
Gas turbine outlet pressure, P; 
height of greenhouse, H 
GT cycle efficiency, Nerep 
Pressure drops, AP 
Organic Rankine cycle turbine inlet pressure, 
Porc, Tur 
Organic Rankine cycle turbine inlet temperature, 
Tore, Tur 
Organic Rankine cycle condenser pressure, Porc 
cond 
Greenhouse desired indoor temperature, Tox, in 
Greenhouse HEX (for winter) outlet temperature, 
Ti7 
Overall heat transfer coefficient for heat 


exchangers 


90 
1.5 


Water-ammonia mixture 


0.4 
0.2 
104 
3 
35 


°C 


kPa 
°C 


°C 


W/m’? K 


4.1. Mode! validation 


A modified MATLAB code was used to solve the current model; the 
results were then for each cycle, in comparison to other articles to 
demonstrate the dependability of the code. Table 7 lists the outcomes of the 
modeling of the gas turbine cycle, the CO2 collecting unit, and the ORC 


subsystem [7,8,24,43]. 


Table 7: Model validation for Brayton cycle and Organic Rankine cycle 


Brayton cycle 


Exergy efficiency (%) 


Organic Rankine cycle 


Total power generated (MW) 


Current 


Ref [59] MAPE (%) 


Modeling 


Model 


Ref [61] 


MAPE (%) 


33.87 


34 


0.38 


51.22 


51.35 


0.25 


Net power production of Gas turbine 


Total cost rate of the system ($ h"') 


cycle (MW) 
Current Reference Reference 
. MAPE (%) Model MAPE (%) 
Modeling [10] [61] 
37.18 37.51 0.87 2.91 2.88 1.04 


A strong coherency between the references and the current model is 
seen, as shown in Table 7. It denotes that the programming code has 
undergone acceptable validation and is now ready to be applied to other 
computations. The exergy efficiency, net power production, and overall 
price of the GT and organic Rankine cycles are therefore smaller than 
typical engineering error, as shown in the relative errors of Table 7. As a 
result, the generated MATLAB programming is validated according to 


published research and is suitable for the current proposed model. 
4.2. Parametric analysis 


On the main goal, such as energy and exergy efficiencies, annual net 
interest and the amount of CO2 released to the ambient, the impact of the 
highest crucial factors, such as greenhouse width, length and height, 
equivalence ratio, grid net electricity required, the organic Rankine cycle 
HTF, and the compression ratio in the GTC, is evaluated. Therefore, each 
factor is viewed as changeable during the others are stable, and each 


output's change will be shown. 
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Figure 4: The variation of generator pressure on the objective functions 


Figure 4 illustrates the variation of generator pressure on the 
objective functions. Energy and exergy efficiencies are increasing with a 
low slope. However, the annual net interest is descending rapidly, because 
the cost function of the absorption refrigerant cycle is a function of the 
pressures deeply. The carbon dioxide emission cost is related to the Brayton 
cycle, which is independent of the absorption refrigerant cycle, so the CO2 
emission is constant. The graphs related to summer are orange and red, and 
those related to winter are blue and green, which overlap because the 
critical state of cold and hot seasons is considered for the calculations of the 
respective system. 
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Figure 5: The variation of greenhouse length on the objective functions 


The influence of greenhouse length on the model functions was 
assessed and is depicted in Figure 5 since the construction of the 
greenhouse and its dimensions are among the primary goals of the current 
article. Nevertheless, the yearly total interest and carbon dioxide 
emission rate have mean amounts each year. Energy and exergy have 
distinct amounts in the summer and the winter. As a result, the first y-axis 
on the left hand of the graph has an orange hue to demonstrate energy 
efficiency, and the graphs are blue and orange in color. Energy efficiency is 
represented by the two green and red graphs, and its label sign is shown in 
red on the right y-axis of the left side. The caption for the purple graph, 
which represents the system's yearly net interest, may be found on the first 
y-axis of the figure's right hand side. Finally, the carbon dioxide cost rate is 
represented in the graphic for the right y-axis of the right by the brown 
color. 


These four purposes and the length of the greenhouse have an 
obvious link. The more a greenhouse's dimensions grow, the more energy 
and exergy are required for it. As a result, as the length of the greenhouse 
increases, so do the energy and exergy efficiencies. Since greenhouse 
production is taken into account while calculating yearly income since more 
items would be produced yearly with a larger greenhouse, the annual 
interest would be more than the costs and the annual net interest would 
rise as the greenhouse's length increased. However, because the 
greenhouse is longer, the rate of carbon dioxide capture would be higher. 
As a result, as shown by the brown hue, the rate of carbon dioxide release 


into the atmosphere will drop and the consequent penalty will reduce. 


1le6 
— 51.6 
=—®— Winter energy efficiency 
—m Summer energy efficiency 51.4 
36 —* Winter exergy efficiency 
—+ Summer exergy efficiency 
51.2 
se se Fe 
XS XS = 
> 35 Es an 51.0 
U U = 
Los = co) 
wv ao) = 
= a 7 50.8 
5 ia © 
> 34 > 2 
{@)) ie) © 
x c r= 
uu WwW <q 
50.4 
33 
50.2 
32 50.0 


50 100 150 200 250 300 350 400 
Greenhouse Width [m] 


Figure 6: The variation of greenhouse width on the objective functions 


The impact of a greenhouse's width on proposed system outputs is 


Total emission costs [$/hour] 


equivalent to that of a greenhouse's length, although the width's range 
differs somewhat from that of the length's range. The maximum amount for 
yearly total interest and the lowest amount for carbon dioxide release 
penalty rate will change since the highest width is 390 m and the less 
length is 500 m. The energy and exergy efficiencies, however, are 
comparable. Figure 6 displays the energy efficiency in summer and winter 
using red and green color schemes, respectively. Energy efficiency charts 
for the summer and winter are illustrated in the figure. Annual net interest 
is shown in a purple diagram, while the carbon dioxide emission cost rate is 
shown in a brown graphic. The justifications and considerations around 
modifications of these criteria for changing width are the same as those for 


changing length. 
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Figure 7: The variation of net gas turbine power on the objective functions 


Figure 7 illustrates how the system's primary outputs change when 
the net power of the gas turbine cycle changes. The total power impacts the 
fuel inflow of the turbine, which is a determining factor for energy and 
exergy efficiency. It is the GT output electricity less’ the 
compressors' required work. The fuel inflow in the combustion chamber is 
modest due to the low power requirements from the cycle of gas turbines. 
Therefore, there is high energy and exergy efficiency. The mass flow rate of 
fuel, which is a result of the total power needed from the GTC, is the key 
word in this study. In addition, as power raises, the cost of carbon dioxide 
emissions ascends and the yearly interest rate increases as more electricity 


is produced and sold. Energy is represented by the blue and orange 


300 


200 


Total emission costs [$/hour] 


diagrams, while exergy is represented by the red and green portraits. The 
yearly net interest is represented by the purple color, and the carbon 


dioxide emission cost rate is shown by the brown color portrait. 
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Figure 8: The variation of equivalence ratio on the objective functions and system 


Figure 8 shows how the primary outputs change when the 
equivalence ratio is changed. The difference between the actual fuel to air 
ratio and the stoichiometric fuel to air ratio is known as the equivalence 


ratio. Because the stoichiometric and actual fuel moles in this relationship 
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are comparable, the relationship will be reduced to the proportion of 
stoichiometric air to real air. In addition, for further assessment and 
analysis, Figure 10 and Figure 11 illustrate the variation of high and low 
pressures in the ORC, which affect the performance of the exchangers, 
pump, turbine, objective functions, and total efficiencies. Figure 10 
illustrates the high pressure of ORC on the objective functions, including 
exergy efficiency, energy efficiency, annual net interest, and CO2 emission 
cost. Exergy efficiencies and energy efficiencies are decreasing for both hot 
and cold seasons by ascending the high pressure of the Organic Rankine 
cycle. The carbon dioxide emission cost is related to the Brayton cycle, 
which is independent of the organic Rankine cycle, so the CO2 emission is 
constant. The annual net interest is decreasing by increasing the pressure 
ratio; this is why the cost of the organic Rankine cycle is soaring when the 
respective pressure ratio is ascending. Figure 11 depicts the low pressure 
of the organic Rankine cycle on the objective functions, including exergy 
efficiency, energy efficiency, annual net interest, and CO2 emission cost. 
The carbon dioxide emission cost is related to the Brayton cycle, which is 
independent of the organic Rankine cycle, so the CO2 emission is constant, 
similar to the high pressure of ORC. Similar to the variation of the high 
pressure of ORC, the annual net interest is decreasing by increasing the low 
pressure; this is why the cost of the organic Rankine cycle is soaring when 


the respective pressure ratio is ascending. 
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Figure 9: The variation of pressure ratio of compressors on the objective functions and 
system 
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Figure 10: : The variation of high pressure of ORC on the objective functions and system 
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Figure 11: The variation of low pressure of ORC on the objective functions and system 


The compression ratio, along with other factors, has an impact on the 
gas turbine's output work, which is its main function. When the 
compression ratio is at a high level, the output power needed by the gas 
turbine will grow, and when it is low, the output power of the GT will drop. 
Thus, it follows that the compression ratio has a direct impact on energy 
generation. Although, because of this factor's increased demand on 
compressors, the net energy output will gradually grow as the compression 
ratio rises. Equations (16) and (17) show that the energy and exergy 
efficiencies increase gradually, as seen in Figure 9. (17). After a certain 
compression ratio, the efficiency variances tend to disappear entirely. In 
lower ratios, the variations are more noticeable. Additionally, when the 
compression ratio increases, the capital price of these parts increases, as 


shown in table 4 and the cost functions for the compressor and GT. As a 


result, as shown by the purple color chart, the cost rate and yearly cost of 
the assessed system both increase, which in turn reduces the net annual 
interest. The gas turbine inlet temperature is the only characteristic that 
the compression ratio impacts; its effect is similar to that of the 
equivalence-ratio. Therefore, as the brown color diagram has demonstrated, 
neither the incoming fuel nor the output carbon dioxide will vary as a result 


of the compression ratio. 
4.3. Outcomes of machine learning method 


Prior to adopting neural networks, the runtime of the optimization 
process was approximated based on the population size and the maximum 
generation possible with the personal computer system in use. It lasted 
around 4 seconds to complete the calculation, and 51.38 hours are expected 
to pass throughout the optimal solution. However, utilizing the trained 
Neural Network, each calculation takes just a tiny portion of a second to 
perform, and the estimated run period process for the applied optimization 


by the neural network is 12 minutes. 


As was mentioned, a set of input and output data is necessary to train a 
supervised neural network. As a result, a set of data was produced using the 
Minitab program and MATLAB code. The input and output data from the 
MATLAB function which was employed to solve the functions of the main 
system were used for the neural network training. Figure 12 depicts the 
trend of training for all types of data and assesses the training performance. 
Figure 13 provides the learning outcomes for each output function. There is 
a gap between data due to the limitations of the Minitab-generated data and 
the lack of data for this gap. The training procedure for yearly interest is 
accurate using a maximum regression, though, and there is a suitable fit 
between the trained inputs and the outputs for the remainder of the data 


set. 
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Best Validation Performance is 3.8159e-10 at epoch 182 
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Figure 12: The trend of traininig for various data 
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Figure 13: Results of regression for input and output data for the main objective functions 
(energy efficiency, exergy efficiency, annual net cost and CO2 emission cost) 


4.4. Optimization results 


Non-dominated sorting genetic algorithm IT (NSGAII) is also a method 
used for multi-objective optimization. Despite the NSGAII being faster 
based on a fast non-dominated sorting algorithm, sharing, elitism, and 


crowded comparison, it does not solve the main problem, which is time- 


consuming run-time. The best way might be to train a precise neural 
network model because utilizing the model is faster than any optimization 
method because optimization through any method is very time-consuming 
compared to a trained model and a conventional optimization method. This 
is why running the trained mode will be performed in a tiny fraction of a 
second. Nevertheless, after training the model, the optimization is 
affordable through conventional methods such as a multi objective 


optimization genetic algorithm for time, energy, etc. 


A genetic algorithm was employed to order to applying an 
optimization method for the presented functions once the neural network 
had been trained. As seen in Figure 14, the Pareto frontier is provided. 
Other depicts show the three dimensional Pareto frontier from various 
angles. Despite the fact that energy efficiency changes from summer to 
winter, the variation is too small to be classified as two separate goal 
functions. The exergy efficiency, yearly total interest, and carbon dioxide 
release cost rate are, therefore, the three most crucial objective functions 
and the exergy efficiency values are based on winter. As demonstrated, the 
Linmap approach was used to determine the optimal position for these 


three goals, which is as follows: 


3 
d= {5 U [f —Fideat,) (32) 
j=l 


The acquired quantities are the most optimal feasible values for each 


objective function when the gain value for d, is minimized. 
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Figure 14: Illustrating Pareto frontier for the three main objective functions 


The Genetic Algorithm (GA) is the optimization method that is employed, 
and Table 8 lists all of the specific characteristics and choices that were 
made for this algorithm. As it is known, in multi objective optimization 
genetic algorithms, there is no single optimum point. Instead, the 
objective is to find a set of solutions that represent the Pareto front. To 
being ensure finding the best solution the optimization process repeated 
at least three times with proper parameters, such as population size, 
mutation rate, crossover rate, and function tolerance. To avoid stalling at 
a not-optimum area and searching enough vicinities of each range, the 
population size was adjusted to 200, and the mutation rate was 
considered to not be very high value. In addition, the major factor to 
assess the convergence criteria and ensuring finding the optimum Pareto 
front is the tolerance of function which is adjusted to 2 «10°. 


Table 8: The optimization adjustments for the Genetic Algorithm 


Multi objective Genetic 


Parameters Algorithm 
Population size 200 
Number of decision 
variables : 
Function tolerance 210 = 
Crossover fraction 80% 
Generations 185 
Max stall generations 65 
Run-time 185 x250 x4|second)=51.38 hou 


As was previously mentioned, a deep neural network was deployed because 
the execution time of pure GA and MATLAB functions was insufficient. 
When utilizing a neural network, the entire optimization process takes 12 
minutes to complete, which is much less time than when utilizing a MATLAB 
function on the same system. Based on Table 9, the model's limitations are 
taken into account. While analyzing the proposed system, various 
constraints, including linear and non-linear constraints were applied to the 
generated codes based on real operating conditions and thermophysical 
limitations to ensure the analysis was correct or not. At first, the generated 
code will be run to check the logical and physical conditions; if so, the 
analysis will continue. By the way, there are several constraint functions in 
the codes to satisfy the real conditions. Table 10 also includes the quantity 
of inputs, the relevant ranges, and the justification for every solution 
variable restriction. 


Table 9:Tthe constraints and limitations of the proposed system 


Parameters Constraints 


Greater than 
70m 


Greenhouse width Greater than 


Greenhouse length 


50m 


: Greater than 
Temperature of the Gas turbine(°C) 


1100 
Output temperature of combustion 
1200 °c 
chamber 
Enthalpy output of the HE of ARC 695.2 kJ 


P30 > P30 saturated 
ORC pressures 
and P30 > P32 


0.5MPa<P< 
ORC high pressures 
10 MPa 
0.1 MPa<P< 
ORC low pressures 
1 MPa 
Enthalpy output of the HE of 
eee 536.7 kJ 


Organic Rankine cycle 


Table 10: Defined decision variables and the lower bound and upper bound 


Lowe Uppe 


r r 
Input variables 
boun boun 


d d 
Generator pressure (kPa) 500 1000 
Length of greenhouse (m) 70 450 
Width of greenhouse (m) 50 450 
Total power of the GT cycle 48,00 
1,500 


(kW) 
Equivalence ratio of the CC (-) 0.12 0.52 
Total pressure ratio of the 


9 21 
compressors in GT cycle (-) 


The Organic Rankine cycle fluid was also examined, and a number of fluids 


were taken into account for the system before the HTF with the highest 


efficiency was chosen for the proposed system, as shown in Table 11. 


Table 11: Various of the Organic Rankine cycle HTF and GT power output 


Organic Rankine Organic Rankine cycle 
cycle working turbine power generation 
fluid (kW) 
Carbon dioxide 1452 
Propane 785 
R134a 547 
IsoButane 751 
Butene 701 
R227ea 789 
R245fa 315 


Table 10 makes clear how important the working fluid is for the power 
output of the Organic Rankine cycle turbine. Table 10 shows various 
Organic Rankine cycle fluids that operate under steady conditions. The ORC 
cycle has a lower pressure of 900 (kPa) and a high pressure of 10000 (kPa), 
respectively. The outcomes show that carbon dioxide outperforms the other 
Organic Rankine cycle fluids that were tested, as shown in Table 11. The 
input and output data used to get these findings are based on the optimum 
results of optimization. Table 12 and Table 13 provide the ideal inputs and 
outcomes, which have been examined to the outputs of the current system 


without the process of optimization. 


Table 12: Factors that affect system decisions for both the old system and the improved 


system 
j . Previous 
Variables optimum value . 
conditions 
P 8792 1000 


generator pressure(kPa) 


Lex (m) 495.24 100 


Wicu (m) 384.63 100 


Wer, net (KW) 8525.17 2000 
@p (-) 0.44277 0.2 
RP omp (-) 12 10 


Table 13: Major outputs of the system prior to and following optimization, together with 
their improved percentage 


After Before 
: ss . oe _ Improvem 
Main parameters optimizati optimizati 
ent (%) 
on on 
Power of the Organic 
2251.2 969 132 
Rankine cycle (kW) 
Total Energy efficienc 
ae - 49.2 45.3 8.6 
(%) 
Total Exergy efficienc 
eh 2 36.71 31.6 16 
(%) 
Total Annual net interest 18,498,671. 
§) 22,522,356 . 21.8 


Table 14 contains a list of the system's defined states based on the optimal 
points. When analyzing and designing the system, it might be useful to 
know the various properties of utilized sub-systems and components and 
other conditions of the proposed system. 


Table 14: Properties of each point according to the optimized condition 


Mass flow 
stat Temperature Pressure . Enthalpy Entropy (kJ/ 
rate 
es (°C) (kPa) (kJ/kg) kg K) 
(kg/s) 
1 20 101.325 19.86 303.5 6.8403 
2 167.3 351 19.86 452.1 6.8952 


3 20 347.5 19.86 302.9 6.4848 


4 168.7 1203.7 19.86 452.8 6.5425 
5 26 101.325 1.4867 -3 0.0496 
6 1131.2 1203.7 20.35 1531.7 7.8202 
7 571.3 104 20.35 881.3 7.9346 
9 571.3 104 20.35 881.3 7.9346 
10 388.1 104 20.35 683.5 4.4853 
11 214.7 103 20.35 501.4 7.3372 
12 211 102 20 497 7.3481 
13 201 102 1.35 679.2 3.1636 
16 91 101.325 2 230.6 0.6969 
17 27 101.325 2 74.6 0.2245 
20 139.7 1350 5.57 538.9 2.0183 
21 96.5 1350 7.1 372 1.8234 
22 36 1350 2.14 1499.3 5.2109 
23 36 1350 2.14 377 1.5684 
24 -12.7 250 2.14 377 1.6489 
25 -12.7 250 2.14 1457.1 5.8113 
26 46.4 1350 5.57 119.2 0.8646 
27 46.6 250 5.57 119.2 0.8684 
28 34.6 250 7.1 55.5 0.8952 
29 34.8 1350 7.1 57.1 0.8962 
30 135 9000 47.5 282.6 1.2241 
31 49 1000 47.5 241 1.2337 
32 49 1000 47.5 197.8 0.955 
33 55 9000 47.5 206.8 0.9616 
34 26 101.325 16 115.9 0.3693 
35 71 101.325 16 304.1 0.9572 
36 26 101.325 11 115.9 0.3693 
37 71 101.325 11 214.6 0.9571 


Table 15 also provides energy and financial streams for the various utilized 


sub-systems and components. 


Table 15: Cost rate and exergy analysis for components 


Exergy destruction 


Component aw) Cost rate ($/h) 
Compressor I 308.8 2.145 
Intercooler 136.4 0.098 
Compressor II 324.3 2.145 
Combustion chamber 10,119 23 


Gas turbine 1,573.2 21.9 


Heat exchanger I 1,626.7 0.012 


Chiller HTF 1,570 6.44 
Greenhouse HEX 368.6 0.012 
Organic Rankine cycle HEX 458.6 2.16 

Organic Rankine cycle 
: 133 10.9 

turbine 

Organic Rankine cycle 

203.5 0.43 
condenser 

Organic Rankine cycle 

91.7 0.22 


pump 


Figure 15 illustrates the depiction of a flow diagram for the rate of 
produced and consumed carbon dioxide by the greenhouse. It is obvious 
that refraining from releasing carbon dioxide into the environment will 
reduce environmental impacts and produce more greenhouse crops. Thus, 
an increase in profit costs is triggered. Flow diagrams, shown in Figure 16 
and Figure 17, depict, respectively, the dissipated energy streams that go 
through the system and the cost per unit of each item. The combustion 
chamber is where the majority of energy is destroyed, as seen in Figure 16, 
which is expected given the elementally of combustion and _ its 
irreversibility. The cost rate of gasoline use in Figure 17 is significantly 
higher than other cost rates. These two numbers might provide an overview 


of the phenomena from an energy and financial standpoint. 


Figure 15: Depiction of flow diagram for rate of produced and consumed carbon dioxide by 
greenhouse 


Compressor]: 308.8kKW ~~ 


Gas turbine: 1,573.2kW | 


Heat exchanger1: 1,626.7kW g 


ondense 


Inlet exergy from the air: 4.0kW c r 
ORC pump: 


Figure 16: Depiction of flow diagram for exergy destruction in the assessed system 


Combustion chamber: 2.300$/hr Compressor]: 2.145$/hr 
Intercooler: 0.098$/hr 
Compressor2: 2.145$/hr 


Fuel: 1,752.000$/hr 
Total cost rate: 1,946.100$/hr cost: 1,946.100$/hr 


Gas turbine: 21.900$/hr 

Heat exchanger1!: 0.012$/hr 
Absorption chiller: 6.440$/hr 
Greenhouse heat exchanger: 0.012$/hr 
ORC Heat exchanger: 2.160$/hr 


ORC condenser: 0.430$/hr 
Other cost rate: 66.338$/hr ORC pump: 0.2203/hr 


Environmental cost rate: 79.000$/hr = ORC turbine: 10.900$/hr 


Figure 17: Depiction of flow diagram for cost rate for each sub-system and elements 


5. Conclusion 


The current article analyzes a real-world case study to suggest, model, and 
optimize a CCHP and greenhouse system using deep learning to improve 
the 4E objectives and carbon dioxide capture. This study also proposes and 
analyzes building or merging a greenhouse near common Gas turbine 
cycles. This research did not focus on separation techniques and catalytic 
converter, and neglected the required energy to do this, but as mentioned 
in section 2 about the separation techniques, the adsorption process is 
preferred for this system because of lower energy consumption and 
financial costs. The process of CO2 separation involves several steps, but if 
a greenhouse is built or merged near the common Brayton cycle, as this 
study proposed, it causes a decrease in the cost and energy of compression 
and transformation of the CO2. The main achievement is to reduce carbon 
dioxide emissions from 0.9183 to 0.41 kgCO2/s to eliminate environmental 
impacts and approach zero carbon dioxide production. Exergy and energy 


efficiencies reached their maximums at 36.71% and 49.2%, respectively. 


The net yearly interest raised by 21.8%, to $22.5 million, as a result of 
greater harvests brought on by the greenhouse's perfect CO2 concentration 
and efficiency. Additionally, the optimization and analysis times were 
decreased substantially by the proposed machine learning method, and an 
accurate prediction was provided. It also includes an air conditioning supply 
component to utilize the power plant's heat loss for the greenhouse and an 
organic Rankine cycle to produce additional electricity load using the power 
plant's energy loss. The following outcomes and implications make the 


current research significant: 


vy A comprehensive energy, exergy, economic, and environmental 
analysis is applied to the proposed system, which shows that this 
approach is substantially adaptable to the other power plants and 
greenhouses in order to approach the elimination of carbon dioxide 
production. 

¥ The primary goal is to lessen the negative environmental effects of 
traditional energy sources, and it has been illustrated that the 
greenhouse absorbed carbon dioxide emissions by around 50%, which 
leads to the refrain of the production of 16,029 tons of carbon dioxide 
annually. 

v¥ By supplying the carbon dioxide of power plants need to grow, selling 
the power produced, and lowering the carbon dioxide emission 
penalty, greenhouse goods have improved, resulting in greater money 
and lower system costs. 

v Also, it is very commendable that a neural network was applied to 
compute the thermodynamic’ properties, economics, and 
the objectives to analysis and optimization in order to shorten the 
computation time and cost, and make accurate predictions. Thus, the 
optimization time improved 257 times after using the trained machine 
learning model. 


¥Y In addition to benefiting greenhouse goods and_ cleaning 


the environment, using greenhouses to absorb carbon dioxide also 


generates economic benefits due to the release of penalties. 


¥ Additionally, using an optimization technique avoids’ certain 
unnecessary costs from arising during the system's development, and 
using the ideal input variables results for the systems in order to work 
at optimal performance and conditions. 
Nomenclature 
A Area of surface (m7) Subscripts 
Acn Cree et Orne 0 Ambient condition 
greenhouse 
B Biases ph physical 
Cc Unit cost ($/kg) Refr Reference 
C, Specific heat (kJ/kg K) SS Strong solution 
C Cost rate ($/s) Tor Turbine 
d Distance of objective function ws Waal solution 
values from reference point 
ex Specific exergy (kJ/kg) Greek Symbols 
EMI1 Extra emission rate of CO, E Heat exchanger effectiveness 
Ex Exergy rate (kW) Q Equivalence ratio 
Ex, | Exergy destruction rate (kW) Deg Maintenance factor 
f,, | Value of ith objective function n Efficiency (%) 
F, Partition function p Energy function 
h Specific enthalpy (kJ/kg) Abbreviation 
H Enthalpy (kJ) ARF Absorption refrigerant fluid 
I Interest rate wt weight 
L Greenhouse length (m) db dry basis 
i Visible layer HTF Heat transfer fluid 
NSGA Non-dominated sorting 
Dy Pedenrayr I genetic algorithm II 


m Mass flow rate (kg/s) NOL Number of needed lamps 
n Life time of system CRF Capital recovery factor 
N Operating years GA Genetic algorithm 
P Pressure (kPa, MPa) HEX Heat exchanger 
QO Heat transfer rate (kW) ORC Organic Rankine Cycle 
R Required CO>2 rate fora CCH Combined Cooling, Heating, 
greenhouse per unit of area and Power 
Ss Specific entropy (kJ/kg.K) CT Cost index 
T Temperature (°C) a Deep Belief Neural Network 
WwW Weight of units GC Gas cycle 
Wi Greenhouse width GH Greenhouse 
W Work production rate (kW) GT Gas turbine 
X Mass fraction HHV Low Heating Value 
Zz Cost ($) LHV Low Heating Value 
AH Change in enthalpy (kJ) ARC | Absorption Refrigerant Cycle 
AS Change in entropy (kJ/K) COP Coefficient of Perform 
AT temperature difference (°C) RBM ae d Boltzmann 
achine 
5 Catia ee Mean sage eee 
Abs Absorber TSA pepe swing 
adsorption 
ch Chemical 
CC Combustion chamber 
a Compressor 
Cond Condenser 
e Exit 
Gen Generator 
Int Intercooler 
Out Outlet 
ph physical 
Refr Reference 
bo yo! Strong solution 
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